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Abstract 

The aim of this paper is to get a complete list of positive definite 
symmetric matrices with integer entries 



solutions to the system of equations 



b d 



such that all complex 



1 X — X-yX^ 
1 - X2 = x\x% 

are real. This result is related to Nahm's conjecture in rank 2 case. 

1 Introduction 

The investigation of torus partition functions of certain conformal field the- 
ories is one of the central topics relating physics and number theory. From 
physical considerations, in general, one expects those functions to have nice 
modular transformation properties. }Zh96| is a fundamental paper devoted 
to establish such expectations rigorously by using the theory of vertex op- 
erator algebras. In [N07| . Nahm considered certain rational conformal field 
theories with integrable perturbations, and argued that their partition func- 
tions have a canonical sum representation in terms of g-hypergeometric se- 
ries. Consequently, he conjectured a partial answer to the question of when 
a particular q-hypergeometric series is modular. This is called Nahm's con- 
jecture, and is discussed in [Z07j . (In particular, one may consult pages 40 
and 41 of this paper for a precise statement of the conjecture.) 
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To state Nahm's conjecture, we consider the g-hypergeometric series 



^n t An+B t n+C 



fA,B,c(z) 



(?)m •••(?) 



n=(m,...,n r )e(Z> ) r 



where (q) n denotes the product (1 — q)(l — q 2 )...(l — q n ), A is a positive 
definite symmetric r x r matrix, B is a vector of length r, and C is a scalar, 
all three with rational coefficients. We will call r the rank throughout the 
paper. Associated with A = (a^) we consider the system of r equations of 
r variables x\, x r 



The definition of (pQ) needs to be made more precise when there are nonin- 
tegral entries in A. However, in this paper we only consider the case when 
all entries are integral, so there is no problem. There are only finitely many 
solutions to ([!]), so all solutions lie in Q. For any solution x = (xi, ...,x r ), we 
consider the element £ a = [x\] + ... + [x r ] G %[F}, where F is the number field 
Q(xi, ...,x r ). £ x defines an element in the Bloch group B(F) by a standard 
construction |Z07| . Then Nahm's conjecture asserts that the following are 
equivalent: 

(i) The element £ x is a torsion element of B(F) for every solution x of (pQ) 

(ii) There exist B and C such that fA,B,c( z ) is a modular function. 

It is not hard to show that Nahm's conjecture holds in rank 1 case (i.e., 
r = 1): one may see |Z07| for example. In particular, A = 1,2 are the 
only integral values of A that satisfy (i). However, from rank 2 and above, 
both directions are open. Obviously, it will be very useful if one can get 
a complete list of matrices A such that condition (i) above holds (Let us 
denote this list by L). Presumably, however, this is a hard task. Instead, 
one may ask an easier question: to determine the set of matrices A such 
that all solutions to ([T]) are real. For any number field F, it is well known 
that the free part of B(F) is isomorphic to 27 2 , where 2r2 is the number 
of complex embeddings of F. So if all solutions to ([T]) are real, then F is 
totally real, = 0, and any element of B(F) is torsion. Therefore for each 
rank r, the set of these matrices is a subset of L. We wish to investigate 
this subset for two reasons: first, it is reasonable to expect that this subset 
constitutes a substantial part of L (From |Z07] and [N07] . one may look at 



r 



l-x i = Y[x a j 



(i = l,...,r). 



(1) 
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the available examples of Nahm's conjecture to get an idea of this. Except 
some trivial infinite families, most examples are in this subset.); second, we 
hope this subset is much more tractable. 

In this paper we focus our attention to rank 2 case, and consider only 
the case when all entries of A are integers. With these restrictions, we will 
identify this subset exactly by using Bezout's theorem. The idea is very 
simple: when r = 2, we introduce a new variable z and consider homoge- 
neous equations corresponding to ([1]). Then Bezout's theorem tells us the 
exact number of complex solutions to the system of homogeneous equations, 
counting multiplicity. There are solutions to the homogeneous system of 
equations, corresponding to z = 0, which are not solutions to ([T]). We can 
estimate the multiplicities of these solutions, and this gives us a lower bound 
for the number of solutions to (PQ), counting multiplicity. The technique we 
will employ to estimate the multiplicities is the method of local analytic 
parametrization. For reference, one can see [W50j . chapter IV, sections 1-5. 
On the other hand, we prove an absolute upper bound for the number of real 
solutions to ([1]), counting multiplicity. Combining these, together with the 
condition that all solutions to (pQ) are real, we obtain inequalities for a, b, d, 
which are sharp enough to enable us to determine all possibilities. 

Our aim is to prove 



Theorem 1.1. If 



a b 
b d 



is a positive definite symmetric matrix with integer 



entries such that all complex solutions to the system of equations 



X2 



x l x 2 



(2) 



are real, then 



a b 
b d 



equals one of 



-1 1 



-1 



-1 



2 
2 



2 1 
1 1 
1 0' 
2 



1 1 
1 2 

2 0' 

1 



or 



4 2 
2 2 
1 
1 



2 2 
2 4 



In section 2, we will assume b > and prove theorem 1.1 in this case. 
After we finish the proof of b > case, it should be very clear how to 
generalize the proof to the case when b is negative (6 = case is trivial 
as it reduces to the rank 1 case). We will include the proof for the latter 
case in section 3. At the end of this paper, we will also indicate a possible 
generalization of the above method to higher rank cases. 



Remark 1.2. Concerning the list of rank 2 examples in [Z07], except the 
trivial infinite family corresponding to central charge one representations 
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of the Virasoro algebra, there is essentially only one known example from 
integral rank 2 case, where all solutions give rise to torsion, but not all 

solutions are real. Namely, the matrix 



"4 


1" 




"1 


1" 






(or 






1 


1 


1 


4 



2 The case when b > 

Without loss of generality, we assume a > d in this section. Thus a > b > 1. 
We consider the following homogeneous system of equations in CP 2 : 



z a+b - xiz^ 6 " 1 - x\x\ = 
z b+d _ X2Z b+d-l _ x b x d = 



(3) 



By Bezout's theorem, the number of solutions to ([3]), counting multiplicity, 
is equal to (b + a)(b + d). Obviously, complex solutions to ([3]) with z ^ 
are in one-to-one correspondence with complex solutions to fl2J). Moreover, 
( j 3 j) has two solutions with z = 0: [x± : X2 ■ z] = [1 : : 0], and [x\ : X2 : z\ = 
[0:1:0]. Let us denote their multiplicities by i\ and 12, respectively. And 
we denote the determinant of the matrix A by A. We have the following 
lemmas for (b + a)(b + d) — i\ —12- 

Lemma 2.1. Let g denote the greatest common divisor of 2b and d. If A = d 
and b > d, we have i\ < b{b + d) + g. 
Otherwise, i\ = min {b(b + d),d(a + b — 1)}. 

Proof. Take x\ = 1, then ([3]) are reduced to 

z a+b _ z a+b-l _ x b = q ( 4a ) 

z b+d - x 2 z b+d - 1 -xi = (4b) 

We would like to take a local analytic parametrization of (|4a|) at z = X2 = 
0, and compute i\ with that. Let e denote the greatest common divisor of 
a + b — 1 and 6, and write a + 6 — 1 = eiti, 6 = e^2- Take 

z = t u \ (5) 
1 

then we have (^-) eM2 = i" 2 — 1. We use (1 — t U2 ) eu 2 to denote its Taylor 
series, then 

x 2 = V , (6) 
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where w is a 6th root of unity, and oj' is any chosen primitive 26th root of 
unity. 

It is straightforward to check that equations ([5]) and ([6]) give a local 
analytic parametrization of (|4ap at z = x 2 = 0, provided that oj varies 
among oj\, where oj\ is a primitive 6th root of unity, and k = 1,2, ...,e. The 
point is that for each point on the affine curve (|4a|) . there exists a unique 
pair of k and t, such that ([5]) and ([6]) give rise to the coordinates of that 
point. For definiteness, let us choose oj\ = (u/) 2 .(Note that these choices are 
not unique. However, any choice will give rise to the same answer for i\ } of 
course.) 

Substituting ^ and ([6]) into (jib]) , then (jibl) becomes 



i i 



j b+d )u 2 _ ri+(w)" 2 (i-r 2 )^q, _ , (l-^)^q, 

1 oj' 1 1 u> ' • 1 ' 

By the theory of local analytic parametrization, i\ equals the sum over k of 
the degrees of the lowest degree terms in t in ([7]). 

Since a > 1, u\ > U2, (6 + d)u 2 < u± + (6 + d — l)u 2 . 

If A ^ d, then (6 + d)u2 ^ du\. Therefore, for each choice of w, the 
degree of the lowest degree term in © equals min {(6 + d)u<2, du\}. So we 
have i\ = min {6(6 + d), d(a + 6 — 1)}. 

If A = d, then (6 + d)u 2 = du\. Since d = A = ad - b 2 , 6 2 = d(a - 1). 
So a — 1 > 6 > d. There are two possibilities: 

(i) 6 = d, then a - 1 = 6. We have % x = I{z 2b+1 - z 2b - x b 2 , z 2b - x 2 z 2b - 1 - 
x b ) = I(z 2b+1 - z 2b - x\, z 2b - x 2 z 2b ~ l - (z 2b+1 - z 2b )) = 6(26) = 6(6 + d), as 
the intersection multiplicity equals the product of multiplicities of the point 
on each curve, if the two curves share no common tangent lines at the point. 

(ii) 6 > d. Then a — 1 > 6. Therefore u% + (6 + d — l)u 2 is greater 

than du\ + u 2 , which is the degree of the second lowest degree term of 

i 

f am( ^---; Consequently, the de gree of the lowest degree term in (171) 

equals (6 + d)u 2 if (^r) d / 1, and equals du\ + u 2 if {^r) d = 1- We write 
d = gd' , 26 = gb' . Since g is the greatest common divisor of 26 and d, d' and 
6' are coprime. Therefore, for oj = oj\ and oj\ = {oj') 2 , {^r) d = 1 only if 6' 
divides 2k — 1, which can happen for at most |^ many k. We then get an 
estimate for iy. 

h<(b + d)u 2 (e - y ) + (dux + u 2 )(-^) = b(b + d) + g. (8) 

□ 
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Lemma 2.2. If A = a, then 12 < b(a + b) + d — 1. 

Otherwise, i 2 = min {6(6 + a), a(d + b — 1)}. 

Proof. The proof is very similar to the proof of the above lemma. Everywhere 
one replaces a,e,ui,U2 by d,f,V\,v 2 , respectively. If A 7^ a and d > 1, then 
> i>2, and the same estimate gives 12 = min {6(6 + a), a(d + 6 — 1)}. If 
A = a, then d > 1, and the difference with the previous case is that now we 
have v\ + V2{a + 6 — 1) < av\ + V2, so the terms of degree v\ + V2(a + 6—1) 
survives, and a rude estimate gives 12 < f{v\+V2{a+b— 1)) = b(a+b) + d— 1. 
If d = 1, then a > A, (a + 6)7^ = f 1 + (a + 6 — \)v2 > av\, and 12 = afv\ = 
a(b + d-l). □ 

Having both lemmas, now let us estimate (a+b)(d+b)— i\— %2- Obviously 
A = a and A = d cannot both happen. We have the following result. 

Lemma 2.3. If A = a or A = d < b, (a + b)(d + b) — i\ — 12 > a — d. 

Otherwise, (a + b)(d + b) — i\ — 12 > a. 

Proof. If A ^ a and one ofA^d, A = d = b holds, then 
(a + b)(d + b) -h -i 2 

= (a + b)(d + b) - mm{b(b + d),d(a + b- 1)} - min {b(b + a), a(d + b - 1)} . 

(9) 

Therefore, 

(a + b)(d + b) -h -ii = A + max{0,(i - A} + max {0, a - A} > a. (10) 

If A = a, then (a + b)(d + b) - i x - i 2 > (a + b){d + b) - b(b + d) - (6(6 + 
a) + (d- 1)) = a-d+ 1. 

If A = d < 6, then (a + b)(d + b) - i x - i 2 > (a + 6)(d + 6) - (6(6 + d) + g) - 
a(d + 6— 1) = a — g > a — d. □ 

Next we prove an absolute upper bound for the number of real solutions 
to ([2]) (counting multiplicity). 

Lemma 2.4. TTie number of real solutions to ([2]), counting multiplicity, is 
at most 9. 

Proof. First of all, there is exactly one solution in (0, l) 2 with multiplicity 
one. (In [Z07], Zagier already mentioned that there is exactly one solution 
in this domain, and moreover this holds in much more general case.) By 
eliminating X2, we get 

(i^i)i + ( i^-l = 0. 

X l 
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In the domain (0, l) 2 , there is no ambiguity on the definition of (jlip . and 
obviously solutions to (|2|) in the domain are in one-to-one correspondence 
with solutions to (jlip in (0, 1). Moreover, for any solution to (jlip given by 
x\ = t, we have the well-defined multiplicity as the valuation of the left 
hand side of (jlip at the point < X\ — r > on the affine Aj.. In other words, 
expand the left hand side of (jlip as a formal power series in [x\ — r), the 
multiplicity of the solution x\ = r equals the degree of the nonzero lowest 
degree term. Furthermore, this multiplicity is the same as the multiplicity 
of the corresponding solution to ([2]) . 

However, since the left hand side of (jlip is strictly decreasing in (0,1), 
one easily sees that there is exactly one solution in (0, 1). Since the derivative 
must be negative, the multiplicity is one. 

Next, we consider solutions outside (0, l) 2 . 

If b is odd, then (jlip is defined unambiguously if we concern only real 
solutions, since there is only one branch of the function x — > xb which maps 
real numbers to real numbers. 

If we have at least one solution with x\ > 1 , then d must be even, and 
X2 < 0. This is case (I). We will show that there are at most 2 solutions in 
this case, counting multiplicity. 

If there exists at least one solution with x\ < 0, then there are two 
possibilities: 

X2 > 0. In this case we have X2 > 1, and a must be even. Exchanging 
indices 1 and 2, we see that from case (I) there can be at most two solutions 
for X2, so at most two solutions for x\ as well. 

X2 < 0. In this case we must have both a and d to be odd. This is case 
(II). We will show that there are at most 6 solutions in this case, counting 
multiplicity. 

Note that any solution to (|2|) with x\ G (0, 1) has to be the unique 
solution in (0, l) 2 , so the above exhausted all possibilities of real solutions. 

If b is even, then we have three possibilities: 
x\ < 0, and X2 < 0. We call this case (III). In this case both a and d 
must be even, and we will show that there are at most 6 solutions, counting 
multiplicity. 

x\ < 0, and < X2 < 1. So a has to be even. We call this case (IV), and we 
will show that there is at most 1 solution, counting multiplicity. 
X2 < 0, and < x\ < 1. So d has to be even. This is the same as case (IV) 
with x\ and X2 switched, so there is at most 1 solution, counting multiplicity. 
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For case (I), (fTT|) can be rewritten as 



,Xl —1.1 ,X\ — 1 . d , 

= -(^H* + "1 = 0. (12) 

X l „, d 

It is easy to see that we need to have A < d, in order that this equation 
has a solution for x\ > 1. We have the following identity for the derivative 
function: 

b{xi - l)xxf'{xx) = x\s d [{d - A)xi + A] - s[a - (a - l) Xl ], (13) 

where 

s = (— ( 14 ) 

(Note that multiplying by invertible elements such as xx — 1, or xi in local 
rings do not affect the multiplicities of solutions. And we are making use 
of the fact that if a polynomial equation has N real solutions in an open 
interval, counting multiplicity, then its derivative has N — 1 real solutions in 
the same interval, counting multiplicity.) 

So f'(xi) can possibly have a solution only when x\ < -Ar. In this case, 
when < xi < f'( Xl ) = iff 



s 



d - 1 x b 1 [(d-A)x 1 + A] =a-(a-l)x 1 . (15) 



But it is easy to see that the left hand side of (|15p is strictly increasing, 
and the right hand side is strictly decreasing. So f'(xi) = has at most 
1 solution counting multiplicity. Consequently, f(x) = has at most 2 
solutions, counting multiplicity. 

For case (IV), by eliminating X2 we have 

,1+X.l,l + X,d 

(— + (^-)*-l = 0, (16) 

X X d 

where x\ = — x, x > 1. 

If A < d, then the left hand side of (|16f) is always greater than 0, so we 
don't have solutions. If A > d, then the left hand side of (fT6|) is strictly 
decreasing, so we have at most one real solution counting multiplicity. 

For case (II) and (III), denote xi = — x, with x > 0, by eliminating xi 
we have 

, . , 1 ~f~ X . 1 , \ ~\~ X ^ d t 

9(*) = -(—^)~ b + (^) ¥ " 1 = °- ( 17 ) 

X x d 
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Table 1: upper bound for the number of real solutions # (for a, b, d, 1 denotes 
odd, and denotes even) 



b 


1 


1 


1 


1 














a 


1 








1 


1 


1 








d 








1 


1 


1 





1 





#< 


1+2 


1+2+2 


1+2 


1+6 


1 


1+1 


1+1 


1+1+1+6 



For the derivative function, we have 

6(1 + x)tg'(x) = -x b sf[(A - d)x + A] + Sl [a + (a - l)x], (18) 

where 

s i = (— ( 19 ) 

So real solutions to g'(x) = are equivalent to real solutions to (counting 
multiplicity) 

(A - d) [ -—^ + d [ -—^- - (a + (a - l)x) = 0, (20) 

where u = v = ^f 2 -. But we have 

(- — ~^~ v - )" = - — ~——^2 — x (quadratic polynomial of x). 

So the second order derivative of the left hand side of ([20]) equals iff 
(1 _|_ x ) u ~ 2 

v= x (a polynomial of x of degree at most 3) = 0. 



So the second order derivative of the left hand side of (j20|) equals has 
at most 3 solutions, g'(x) = has at most 5 solutions, and g{x) = has 
at most 6 solutions, all counting multiplicity. Combining all the above, and 
enumerate all 8 cases of a, b, d being even or odd as in table 1, the lemma is 
proved. □ 

Now let us prove theorem 1.1 for the case b > 0. 

Proof. In [Z07J, they first searched for matrices A with integral entries (ac- 
tually they did it for matrices with rational entries with certain bounds on 
numerator and denominator) whose absolute values are less than or equal 
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to 100 such that L(£) G vr 2 Q, where £ = [x±] + [x2\ is the Bloch group 
element corresponding to the unique solution x\,X2 to ([2]) in (0, l) 2 , and 
L(£) = L(x\)+L(x2), where L[x) is the Rogers dilogarithm function. Among 
these they then identified all matrices such that condition (i) of Nahm's con- 
jecture holds. On the other hand, as we have explained in the introduction, 
for any matrix A satisfying the condition of theorem 1.1, A satisfies con- 
dition (i) of Nahm's conjecture, and in particular G vr 2 Q, which fol- 
lows from the well-definedness of the regulator map (for details one may see 
[Z07J). Meanwhile, combining lemmas 2.3 and 2.4, we have a < 9, or A = a, 
a — b<9,orA = d<b, a — b<9. From each of these conditions one easily 
derives that a,b,d < 100, so A has to be in the Zagier's list. Therefore, only 



the matrices 



"2 f 




and 


1 1 



survive. 



□ 



3 The case when b is negative 

We will prove the following 



Proposition 3.1. if 



a b 
b d 



is a positive definite symmetric matrix with 



integer entries such that all complex solutions to the system of equations 



'1*2 



X*2 — XyX% 



are real, and b is negative, then a, —b, d < 20. 

This proposition will be a consequence of the following two lemmas 3.2 
and 3.3, and it immediately implies the full theorem 1.1, again as there are 
only finitely many remaining cases to check to identify all matrices among 
these such that all solutions to ((2J) are real, and all these cases are already 
checked in [Z07j . The proof will be brief since it is similar to the proof of 
the case b > above. 

Let us write b = — c, with c > 0. Without loss of generality, we assume 
a > d. Then a > c + 1. We consider the system of equations 



x\{\ 



x?(l 



X2) 



x'o 



(21) 



Except for x\ = X2 = 0, solutions to ([2T]) are in one-to-one correspondence 
with solutions to ([2]). Using the method of local analytic parametrization, 
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one easily sees that the multiplicity of this solution is c 2 : write a = h±a±, 
c = h\Ci, where hi is the greatest common divisor of a and c. Then x% = t ci , 

x 2 = — * 10 i , where 9 is a primitive cth root of unity, and k varies 

(i-t c i) fc i c i 

among 1, 2, hi defines a local analytic parametrization of x\ — x 2 (l — x\) 
at x\ = X2 = 0. Substituting this parametrization into x\{\ — X2) = x d while 
keeping in mind that ad > c 2 , one sees easily that for every k, the degree of 
the lowest degree term is c\c. So the multiplicity of this solution is c 2 . 

Lemma 3.2. The number of solutions to (12 1 [) other than x\ = x 2 = 0, 
counting multiplicity, is at least a — 1. 

Proof. Let us denote the quantity in the above lemma by n. We introduce a 
new variable z, and consider homogeneous equations corresponding to (12ip . 
Since a > c + 1, we need to discuss two cases: 
Case (i): d > c + 1 

In this case, the homogeneous equations are 

x%(z a ~ c - x x z a - c - 1 ) =x\ 

A A 1 \ ( 22 ) 

x\{z d - c - x 2 z d - c - 1 ) =x\ 

If a > c + 1, then z = implies x\ = x 2 = 0. So we don't have solutions 
to (I22p with z = 0. Thus, Bezout's theorem implies that n > ad — c 2 = A > 
a -I. 

If a = c + 1, then d = c+ 1, (|22l) becomes 

zg(ai-xi) = x' i ; +1 
- x 2 ) = ^2 +1 

There are c distinct solutions with z = 0: [xi : X2 : ^] = [1 : : 0], 

7ri-\-2kTri 

where Wk = e <= ,k = 0, l,...,c — 1. We calculate the local valuation 
of the function x\z — x\ — 1 on the curve z — x 2 — x^ 1 = 0. In fact we 
have x\z — x 2 — 1 = — {x 2 + l)(x2 + 1). So V{x 2 z — x 2 — 1) equals 2 if 
Wk = —1, and 1 otherwise. There is at most one k such that = — 1, 
so the sum of multiplicities of these solutions is at most c + 1. Therefore 
n > (c + l) 2 - (c + 1) - c 2 = c = a - 1. 
Case (ii): d < c 

In this case, the homogeneous equations are 



x c 2 {z a ~ c - x x z a ~ c - 1 ) = x\ 



X2 ) = x d 2 z c+1 - d ^ 
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We have only one solution with z = 0, namely [x% : xi : z] = [0 : 1 : 0]. Let 
us denote the multiplicity of this solution by i. 

Again we use the method of local analytic parametrization. Let X2 = 1, 
([23]) are reduced to 

z a ~ c - x 1 z a ~ c - 1 - x? = (24a) 
x\(z - 1) - z c+1 - d = (24b) 

Write c = hiC2, c + 1 — d = h^di, where hi is the greatest common divisor 
of c and c + 1 — d. Then z = t C2 , x\ = — t d ' 2 9 k ; where fc varies among 

(t c 2-l) ft 2 c 2 

1,2, ...,/i2 defines a local analytic parametrization of (|24b|) at xi = X2 = 0. 
Substituting this parametrization into (|24ap . we get 

tC 2 {a - c) _i _j_ v ^ (25) 

UC2 _ 1) h 2 c 2 (pCH - 1) ^2=2 

If d > 1, then 02(0 — c) > <i2 + C2(a — c— 1), and we have two possibilities: 

(i) a + d 7^ A + 1, then cfe + C2(a — c — 1) ^ ac?2, and therefore i < c + 1 — 
d + c(a — c — 1) = ac — c 2 — d+1 < c(a — c), 

(ii) a + d = A + 1, then efo + C2(a — c — 1) = acfe > C2(a — c) — C2, so the term 
^c 2 (a-c) surv ives, and i < c(a — c). 

Either case, we have n > a(c + 1) — c 2 — c(a — c) = a. 

\{ d = 1, the second lowest degree term in t d2+C2(a c i 1) 6' fc surv i veS; anc i 

we have i < a — c + 1 < a, and n > a(c + 1) — c 2 — a = c(a — c) > a — 1. □ 



Lemma 3.3. The number of real solutions to (12ip other than x\ = X2 = 0, 
counting multiplicity, is at most 19. 

Proof. Except for the solution in (0, l) 2 , again we have four possible cases, 
identical to the case when b > 0. In each of these cases, the same argument as 
in the cases (II) and (III) for b being positive works, which produces an upper 
bound 6 for the number of real solutions in each case. Again by enumerating 
all 8 cases of a, c, d being even or odd as in table 2, we get an upper bound 
for the number of solutions counting multiplicity to be 1 + 6 + 6 + 6 = 19: 

□ 

Combining lemmas 3.2 and 3.3, lemma 3.1 is thus proved, and our proof 
for theorem 1.1 is complete. 
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Table 2: upper bound for the number of real solutions # (for a, c, d, 1 denotes 
odd, and denotes even) 



c 


1 


1 


1 


1 














a 


1 








1 


1 


1 








d 








1 


1 


1 





1 





#< 


1+6 


1+6+6 


1+6 


1+6 


1 


1+6 


1+6 


1+6+6+6 



Remark 3.4. It looks likely that this method can also be generalized to 
deal with the case when A is a nonintegral matrix. In this case one has to 
formulate the problem a bit more carefully and do all the above analysis 
with more patience. In other words, we expect the above method to provide 
a complete list of rational matrices A such that all solutions to ([2]) are real. 
Lastly, it is conceivable to speculate that the method may also be generalized 
to higher rank cases: Again, our goal is to get some inequalities on matrix 
entries from Bezout's theorem, which hopefully are sharp enough to give an 
upper bound for matrix entries depending only on the rank, thus giving a 
finite set of possibilities for the matrix A for each rank, such that all solutions 
to (pQ) are real. Thanks to the work of A.G. Khovanskii, we already have an 
almost satisfactory substitute for lemma 2.4 in higher rank cases: corollary 7 
on page 80 of the book |K91j . It provides an upper bound for the number of 
nondegenerate real solutions to (pQ) depending only on the rank. (Although 
in rank 2 case, this upper bound is too large compared to ours in lemma 
2.4. But the most important thing is the existence of such an upper bound 
in general.) Consequently, if we can find a good way to estimate the sum of 
multiplicities of solutions to certain homogeneous equations corresponding 
to ([I]), which are not solutions to thus getting a substitute for lemma 
2.3, we may achieve our goal, (provided that we also know how to deal with 
the exceptional cases when some of the real solutions to ([!]) are degenerate) 
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